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We present a measurement of the branching fractions for the Cabibbo-favored radiative decay, 

7, and the Cabibbo-supprossod radiative decay, D'^ (pj. These measurements are 
based on a data sample corresponding to an integrated luminosity of 387.1 fb~^, recorded with 
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the BABAR detector at the PEP-II e e asymmetric-energy coUider operating at center-of-mass 
energies 10.58 and 10.54 GeV. We measure the branching fractions relative to the well-studied 
decay ^ K'-k+ and find B{D° K*°-i)/B{D° K--K+) = (8.43 ± 0.51 ± 0.70) x 10"^ and 
cf>j)/B{D° K--K+) = (7.15 ± 0.78 ± 0.69) x 10"*, where the first error is statistical and 
the second is systematic. This is the first measurement of B{D^ K*^'y). 

PACS numbers: 13.25.Ft, 12.38. Qk, 12.40Vv, ll.30.Hv, 13.20.Fc 



In the 6-quark sector, radiative decay processes have 
provided a rich field in which to study the Standard 
Model of particle physics. Decays such as B ^ pj 
have yielded measurements of the Cabibbo-Kobayashi- 
Maskawa matrix element \Vtd I [J, Q. These decays are 
dominated by short-range electroweak processes, whereas 
long-range contributions are suppressed. The situation is 
reversed in the charm sector, where radiative decays are 
expected to be dominated largely by non-perturbative 
processes, examples of which are shown schematically in 
Fig. [1] Long-range contributions to radiative charm de- 



Mode 


Experimental 


Theoretical [3, 4, 5, 6, 7, 8, i 




B.F. (xlO-^) 


B.F. (xlO"'') 






[10] 0.1-3.4 


D« 


K*°'y < 76 (90% C.L.) [llj 


7-80 


D° 


p°7 < 24 (90% C.L.) [llJ 


0.1-6.3 


D« ^ 


Lo-y < 24 (90% C.L.) [11] 


0.1-0.9 



TABLE I: The current experimental status and theoretical 
predictions for the branching fractions (B.F.) of radiative 
charm decays with vector mesons. 
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FIG. 1: Feynman diagrams for the long-range electromagnetic 
contributions to D° ^ V-y, V = Fi gures (a) and (b) 

show sample vector dominance processes, while (c) and (d) 
are examples of pole diagrams, where the circles signify the 
weak transition and P represents a pseudoscalar meson. 



cays are expected to increase the branching fractions for 
these modes to values of the order of 10~^, whereas short- 
range interactions are predicted to yield rates at the 10^^ 
level. Given the expected dominance of long-range pro- 
cesses, radiative charm decays provide a laboratory in 
which to test these QCD-based calculations. 

Numerous theoretical models have been developed to 
describe these radiative charm decays 
The two most comprehensive studies 
similar amplitudes for the dominant diagrams shown in 
Fig. [TJ The first paper bases predictions on Vector Me- 
son Dominance (VMD) calculations, while the second 
paper uses Heavy-Quark Effective Theory in conjunc- 
tion with Chiral-Lagrangians. Though each approach 
arrives at similar estimates for the magnitudes of the 
individual decay amplitudes, Ref. Q predicts that the 
pole diagrams, shown in Figs. [1] (c) and (d), interfere 
destructively and cancel nearly completely. Ref. 
makes no such predictions. Precise measurements of 
B{D'^ — > Vj, V ~ K*°, cj)) may provide insight into the 



developed to 
predict very 



amount of interference between pole diagrams. 

The first observation of a radiative, but color- 
suppressed, decay process was made by the Belle 
collaboration with a measurement of B{D^ — + (fry) — 
(2.43t°j?(stat.)t°;}^(sys.)) x IQ-^ CLEO II con- 

ducted searches for other radiative decays and established 
the current upper hmit of B{D° K*°j) < 7.6 x 10"'' 
at 90% confidence level (C.L.), as well as upper limits on 
B{D° p°-/) and 6(D° -> ^7) [lH. Table U summarizes 
theoretical predictions and current experimental results. 

In this paper we present the first observation of the 
Cabibbo-favored radiative decay K*^^, as well as 

an improved branching fraction measurement of the pre- 
viously observed decay The analysis is based 
on 387.1 fb~^ of data recorded by the BABAR detector at 
the PEP-II e+e~ asymmetric-energy collider operating 
at center-of-mass (CM) energies of y/s — 10.58 GeV and 
10.54 GeV, and uses approximately 5 x 10^ e+e" cc 
events. 

The BaBAR detector is described in detail else- 
where ■ Charged particle momenta are measured with 
a 5-layer double-sided silicon vertex tracker and a 40- 
layer drift chamber. Charged hadron identification is pro- 
vided by measurements of the specific ionization energy 
loss, dE/dx, in the tracking system and of the Chercnkov 
angle obtained from a ring-imaging Cherenkov detec- 
tor. An electromagnetic calorimeter consisting of 6580 
CsI(Tl) crystals measures shower energy and position for 
electrons and photons. These detector elements are lo- 
cated inside, and coaxial with, the cryostat of a supercon- 
ducting solenoidal magnet, which provides a 1.5 T mag- 
netic field. The instrumented flux return of the magnet 
allows discrimination between muons and pious. 

A detailed Monte Carlo (MC) simulation of the BaBAR 
detector based on GEANT 4 13] is used to validate the 
analysis and determine the reconstruction efficiencies. 
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We optimize our selection criteria using simulated events two photons, one of which is combined with the K 



by maximizing significance, defined as Ns / V Ns + Nb , 
where Ns and Nb denote the number of signal and back- 
ground candidates in the MC simulation. We recon- 
struct radiative — > V^, V = K*^,(l) decays using the 
charged decay modes of the vector meson, K*'^ K^tt^ 
{(j) K-K+) [11]. We form K*^ (0) candidates from 
pairs of oppositely charged tracks identified as K~t:^ 
(K^K^) using the Cherenkov angle measurement of the 
DIRC and dE/dx measurements from the tracking sys- 
tem, and accept any K~tt^ {K~ K^) candidates with 
invariant mass in the range 0.848 to 0.951 GeV/c^ (1.01 
to 1.03GeV/c^). The charged track candidates are fit to 
a common vertex, and a fit probability greater than 0.1% 
is required. 

A photon candidate is defined as an energy deposit in 
the EMC that is not associated with the trajectory of a 
charged track, and which exhibits the expected shower 
shape characteristics. Each such candidate is required to 
have CM energy greater than 0.54 GeV. The charged- 
particle vertex is assumed to be the production point 
of the photon. We suppress the significant background 
from tt" 77 decays by rejecting a photon candidate 
which, when paired with another photon in the event, 
results in an invariant mass consistent with the 7r° mass, 
(0.115 < M(77) < 0.150) GeV/c2. 

Background from random D''^ Vj candidates is re- 
duced by requiring that the candidate be a product of 
the decay D*~^ D^tt"'". A D*~^ candidate is formed by 
combining a candidate with a low-momentum charged 
pion, denoted as irf . These pion candidates are required 
to have CM momentum less than 450MeV/c. We cal- 
culate the mass difference, AM = M{V^Trf) — M{V^) 
and require (0.1435 < AM < 0.1475) GeV/c^. The AM 
distribution of candidates arising from signal decays is 
well-described by a Gaussian distribution function. Our 
selection corresponds to a six-standard deviation interval 
centered on the mean of the Gaussian, and hence retains 
almost all of the signal candidates. We reduce combina- 
toric background from BB events to a negligible level by 
requiring that the CM momentum of the D*'^ candidate 
be greater than 2.62 GeV/c. 

The dominant background in our sample of 
.ftr*°7 candidates results from Z)" X^vr+vr" decays, 
where one of the photons from the tt" decay is paired with 
the kaon and pion from the D'^ decay to closely mimic 
the signal mode. As described above, we use a 7r° veto to 
suppress such events but, given the large branching frac- 
tion of this mode, B{D° K~tt+tt°) = (13.5 ± 0.6)% 
[l5.] , a significant number of such candidates survives. We 
can separate this background from signal on a statistical 
basis because of differences in the K'n^^ invariant mass 
distribution. The background distribution peaks slightly 
below the nominal D'^ mass, and has a different shape 
from that of signal events. An additional background 
arises from —> K*^ri events where the ry decays to 



pair to form an invariant mass within our 13° mass win- 
dow. This contribution peaks well below the nominal 
mass, and it can be separated easily from correctly re- 
constructed — > K*^^ decays. 

The impact of both D° K*°tt° and K*°ri is 

further reduced by using the helicity angle 9h- The 
helicity angle is defined as the angle between the mo- 
mentum of the K*^ meson parent particle {D^) and the 
momentum of the K*^ daughter kaon as measured in the 
K*^ rest frame. Due to angular momentum conservation, 
dN/ dcosOn for signal candidates varies as 1 — cos^ 6h, 
whereas for K*'^tt'^ (ry) events the cosine of the he- 

licity angle is cos^ Oh distributed. The cos 6h distribu- 
tion of 13° — > _ftr^7r+7r° candidates is complicated by the 
interference and overlap of resonant structure in the final 
state Dalitz plot. Based on a MC study an asymmetric 
selection of —0.30 < cosOh < 0.65 is chosen to maximize 
the signal significance. 

Similarly, but to a lesser extent, the signal of the 
Cabibbo-suppressed radiative decay 13° 07 is ob- 
scured by backgrounds from 13° — > 07r° and Z3° (prj 
decays. Due to the small width of the (f> meson, back- 
ground from 13° K^K^n^ transitions with a K^K^ 
invariant mass in the <j) region yields a negligible con- 
tribution to 13° [l6| . Since angular momentum 
conservation dictates that the cosine of the helicity angle 
of the remaining Z3° — > (/)7r° events follow a cos^ 9h dis- 
tribution, we replace the tight cos 6h selection criterion 
used in the 13° K*^^ case with the looser requirement 
jcos^ifl < 0.9 and include cosOh as a variable in the 
fitting procedure. This retains a larger fraction of signal 
events, and so reduces statistical uncertainty. 

We consider other radiative decays which might re- 
flect into the M(^*°7) and M{(j>j) invariant mass dis- 
tributions. Background to Z3° _ftr*°7 may arise from 
13° — > 07 if a kaon from — > K^K^ is mis-identified 
as a pion. Background from 13° p°7 may arise if a 
pion from p° — )■ 7r^7r+ is misidentified as a kaon. Real 
Z3° ^*°7 events can reflect into the M{(f>j) distribu- 
tions if a 7r+ is misidentifled as a . Using MC simula- 
tions, all of these background contributions arc found to 
be negligible. 

We extract the D° -> K*""f yield using an unbinned 
extended maximum likelihood method (E-MLM) to flt 
the M(_R'*°7) invariant mass spectrum. The yield of 
Z3° 07 events is extracted using an E-MLM to flt the 
two dimensional distribution of invariant mass, M(07), 
and helicity, cos 9h ■ 



We use a Crystal Ball (CB) line shape [17[ to model the 
invariant mass distributions for 13° — > i^*°7 (Z3° — + (pj) 
signal events, and background reflections from 13° 
K~7T~^7T^ (13° 07r°) decays. The invariant mass dis- 
tributions of Z3° K*^ri and 13° — > 07y background 
events are modeled with a Gaussian function and a flrst 
order Chebychev polynomial. The remaining combina- 
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toric background decays are modeled with a second order 
Chebychev polynomial. In the <^7, case the cos 0h distri- 
butions of D° — > ipj, — > (/)7r", D'^ — * (prj, and combi- 
natoric background events are all modeled using second 
order Chebychev polynomials. The parameters of these 
probability distribution functions (PDFs) are obtained 
using simulated events and subsequently fixed when fit- 
ting the data. 

We validate the invariant mass PDFs using data. To 
verify that the MC correctly simulates the backgrounds 
and the effects of the missing photon from the tt*^ decay, 
we search our data sample for candidates. 
Since this decay is forbidden by angular momentum con- 
servation, the candidates surviving our selection criteria 
are all combinatoric background or due to 
or — + K^r] decays. 

We select candidates from pairs of oppositely 
charged tracks identified as pions. The pions are re- 
quired to share a common production vertex and have 
an invariant mass in the range (0.490 < Af(7r+7r~) < 
0.505) GeV/c^. Selection criteria for the photon momen- 
tum, £)*+ momentum, AM, and tt*^ veto are identical 
to those used in the V-f analyses. The resulting 

Kgj invariant mass spectrum is fit with a linear combi- 
nation of three PDFs. The first PDF is used to model 
Z?" Kgir'^ candidates, and has the same functional 
form as the one used to model D'^ K^tt^tt^ candi- 
dates. The second PDF is used to model Kgi] can- 
didates, and has the same functional form as that used 
to model D° K*°-q candidates. The third PDF is a 
second order Chebychev polynomial used to model com- 
binatoric background candidates. The shapes for both 
Kgi] and combinatoric background candidates are 
fixed using MC. The Kgir^ signal shape is allowed 
to float in the final fit. Both MC and data are fit in this 
way and we find good agreement. The observed differ- 
ences in the fit parameters are used to correct the CB 
line shape PDFs as described below. 

A second test is performed using D'^ K*'^j can- 
didates taken from the sideband regions defined by 
\ cos 9h\ > 0.9. Very few D° K*°^ candidates are 
seen within this region, leading to a clean sample of 
K~'k^'it'^ decays. The resulting invariant 
mass spectrum is fit using a procedure similar to the 
one used for signal region D° — > K*^^ candidates. The 
only differences are that the K*'^j contribution 

is fixed to zero, and the X^tt+tt" signal shape 

is allowed to float freely. The M {K*°j) distribution of 
D° — > if "TT+TT*^ events is compared between data and 
MC and we find good agreement. 

Potential differences between the D'^ invari- 
ant mass distributions for data and MC are evaluated 



by using 



ifgTr" events. The selection criteria 
for A'^ mesons are identical to those applied in the 
— > analysis. The requirements on AM and 

CM momentum are as before. A 7r° candidate 



consists of a photon pair with invariant mass satisfying 
(0.110 < M(77) < 0.150) GeV/c2, and resultant labora- 
tory momentum greater than 0.540 GeV/c. This resulting 
sample of — > K^n^ candidates is fit to a CB line shape 
and a linear background. 

We used the average difference between the CB line 
shape parameters in MC and these data control samples 
to modify the PDF parameterizations used in the fit. 

The fit results from data and expected signal and back- 
ground contributions from MC are shown in Fig. [2l(a- 
c). The event yields obtained from the E-MLM fit for 
both D° -> 07 and D° ^ K*°-f are N{D° ^ <h\ ^ 



242.6 ± 24.8 and N{D^ 



K* 



'T, 



K* 



K- 

K^TT^) = 2285. 8± 113.2. The reconstruction efficiencies, 
determined using MC, are found to be e(Z?" ^7; (p ^ 



7; 



K* 



K-K+) = (10.8 ± 0.1)% and e{D^ K' 
K-TI+) = (6.4 ±0.1)%. 

In order to avoid uncertainties in the overall normal- 
ization we measure the branching fraction of the ra- 
diative decays relative to B{D^ A'^7r+). We pre- 
pare a K^TT^ datasct following procedures sim- 
ilar to those described above, and find a yield N{D'^ 
K-TT+) = (335.1 ± 4.0) X 10^ with an efhciency of 
6(1)0 ^ (5.3 ± 0.2)%. 

We perform several consistency checks. Our result is 
compared to the cos 9h distribution expected for 
K*^^ by refitting the data in intervals of cos Oh and mea- 
suring N{D^ K*^^) in each interval. The normalized 
and efficiency-corrected result, shown in Fig. 2(d) com- 



pares well to the expected distribution. As an additional 
check, we divide the data into five distinct samples, one 
for each PEP-II run period, and perform the analysis on 
each subset independently. We see a K*'^j signal 

for each run period, and find that the branching ratios 
are consistent within statistical uncertainties. 

We evaluate the systematic uncertainties associated 
with our measurement in several different studies. Sys- 
tematic effects due to the PDF parameterizations of sig- 
nal and backgrounds are determined by generating an 
ensemble of 1,000 random numbers drawn from a normal 
distribution for each PDF parameter, including their cor- 
relations obtained from our fits. We refit the data using 
each of the 1.000 sets of random numbers. The result- 
ing distribution of A^(_D0 Vj) is fit to a Gaussian 
function and the percent standard deviation is taken as 
the systematic error, 5.9% for — > and 4.4% for 

Correcting the D° V-/ and Vt:^ PDF pa- 

rameters using the data control samples induces a sec- 
ond systematic uncertainty in the parameterization of 
the signal shapes. We estimate this effect by indepen- 
dently applying the corrections obtained using each of the 
three control samples. The largest percentage variation 
in N{D'^ — > Vj) is taken as the systematic uncertainty 
associated with this correction; this leads to systematic 
uncertainities of 3.0% for N{D^ — > ^7) and 4.3% for 
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I-. .1 D ^ (i r Background 

BG Fit ^ - 

■ D Tt +BG Fit 

Total Fit 

Data ~ 




1.95 2 

M(<Ixy) (GeV/c=) 



(a) The ^7 invariant mass distribution. 



Generic Background (BG) 
d" -> Ji° Background 
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(b) The 07 hehcity angle distribution. 
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D ^ K T[ Background 
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(c) The K*^'j invariant mass distribution. 



C0S{8„) 

(d) The D° K*°^ hclicity angle 
distribution. 

FIG. 2: Invariant mass and cos 6h distributions for data 
(points) and simulated events (histograms). The curves show 
the fit results and the individual signal and background con- 
tributions. BG refers to the combinatoric background. 



We quantify the difference in particle identifica- 



Systematic cj"(Z^" 


ch-y) (%) a(D" - 




Tracking, vertexing 


1.2 


1.0 


Particle ID 


2.9 


1.1 


7 reconstruction 


1.8 


1.8 


Tv" veto 


1.8 


1.8 


PDF parameter 


5.9 


4.4 


Correcting Vno^v-y 


3.0 


4.3 


and V^Q^v-wO 






Ref. mode efficiency 


1.5 


1.5 


Selection criteria 


5.4 


4.5 


Total systematic effect 


9.6 


8.3 



TABLE II: Summary of all systematic errors for each D" de- 
cay mode. The total systematic uncertainty is obtained by 
adding the individual systematic estimates in quadrature. 



tion (PID) efSciency between data and simulation by 
means of a high-purity control sample of D*~^ —^ D^tt"*", 
D'^ K~iT'^ events, which we divide into intervals of po- 
lar angle and momentum. The change in yield when PID 
selection criteria are applied is computed separately for 
data and for simulated events and the difference is taken 
as a correction factor for that interval. We then weight 
the correction factors according to the expected momen- 
tum and polar-angle distributions of the K*'^j sig- 
nal. While a portion of the PID systematic uncertainty 
for our signal modes is canceled when measuring the 
branching fractions in ratio to D'^ — s- if "tt"*", the residual 
uncertainty is found to be 2.88% for (j>j and 1.10% 
for D" K*°-f. By measuring B{D" -> K*"-f) and 
B{D'^ 07) with respect to 1?° K^ir^ , first-order 
effects from charged particle tracking also cancel, leaving 
only a second order systematic uncertainty of 1.00% for 
£)0 ^ ^*o^ events and 1.20% for cj)-f. We sum- 
marize all systematic uncertainties in Table HH 

In this paper, we report our observation of the 
Cabibbo-favored, but color-suppressed, radiative decay 
K*^^. We also present confirmation of the pre- 
vious measurement of the Cabibbo-suppressed radiative 
decay B{D^ — > (f)"f), but with reduced statistical uncer- 
tainities. The measured branching ratios are 



B{D° -> 07) 
B{DO K-Tr+) 
B{D° K*"-f) 
B{D" K-Tr+) 



= (7.15 ±0.78 ±0.69) x lO"'^ 
= (8.43 ±0.51 ±0.70) X 10"^ 



where the first uncertainty is statistical and the sec- 
ond is systematic. Using the current world average 
of B{D^ ^ K-Tr+) = (3.82 ± 0.07)% [H we obtain the 
following absolute branching fractions: 



6(i?° -> 07) 
6(D" ^ K*°j) 



(2.73 ±0.30 ±0.26) x 10"^ 
(3.22 ±0.20 ±0.27) x 10"*. 
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These results are consistent with the theoretical expec- 
tations of Table HI 

In the context of the vector dominance model the 
largest contribution to radiative decays is expected 
to come from a virtual coupling directly to a sin- 
gle photon, leading to the prediction that the branch- 
ing ratios B{D°_-^ (f)-f)/B{D° K*^) and B{D° 
(l)p°)/B{D° K*°p°) should be equal % Comparing 
our measurements of the radiative D° decays with the 
current world averages [ist we find 



(6.27 ±0.71 ±0.79) x 10"^ 
(6.7 ± 1.6) X 10^2 



in agreement with this prediction. 

If we assume all contributions are from VMD type pro- 
cesses and under the assumption that the p° meson is 
transversely polarized, as has been confirmed experimen- 
tally for D° if*OpO expect B{D'^ V-f) « 
aEMB{D° ~* Vp°) y, where asM = 1/137 is the fine 
structure constant. Using our results we find 



B{D'> ^ 07) 



(0.021 ± 0.005) B{D° 
(0.020 ±0.003) S(L»° 



which in both cases is about a factor of three larger than 
the VMD prediction. This indicates that we are seeing 
enhancements from processes other than VMD, which 
might be explained by incomplete cancellation between 
pole diagrams. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the comput- 
ing organizations that support BaBAR. The collaborat- 
ing institutions wish to thank SLAC for its support and 
kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), CEA and CNRS-IN2P3 
(France), BMBF and DEC (Germany), INEN (Italy), 
EOM (The Netherlands), NFR (Norway), MES (Russia), 
MEC (Spain), and STFC (United Kingdom). Individuals 
have received support from the Marie Curie EIF (Euro- 
pean Union) and the A. P. Sloan Foundation. 



tt 
[1 

[2: 

[s: 
[4: 

[6; 

[7 

[8 

[9: 
[10: 

[11 

[12 
[13 

[14 

[15 

[16 
[17 



Deceased 

Now at Temple University, Philadelphia, Pennsylvania 
19122, USA 

Now at Tel Aviv University, Tel Aviv, 69978, Israel 
Also with Universita di Perugia, Dipartimento di Fisica, 
Perugia, Italy 

Also with Universita di Roma La Sapienza, 1-00185 
Roma, Italy 

Now at University of South Alabama, Mobile, Alabama 
36688, USA 

Also with Universita di Sassari, Sassari, Italy 

A. All and A. Y. Parkhomenko, Eur. Phys. J. C23, 89 

(2002) . 

A. Ali, E. Lunghi, and A. Y. Parkhomenko, Phys. Lett. 
B595, 323 (2004). 

B. Bajc, S. Fajfer, and R. J. Oakes, Phys. Rev. D51, 
2230 (1995). 

B. Bajc, S. Fajfer, and R. J. Oakes, Phys. Rev. D54, 
5883 (1996). 

G. Burdman, E. Golowich, J. L. Hewett, and S. Pakvasa, 
Phys. Rev. D52, 6383 (1995). 

H. -Y. Cheng et al., Phys. Rev. D51, 1199 (1995). 

S. Fajfer, A. Prapotnik, S. Prelovsek, P. Singer, and 
J. Zupan, Nucl. Phys. Proc. Suppl. 115, 93 (2003). 
S. Fajfer and P. Singer, Phys. Rev. D56, 4302 (1997). 
S. Fajfer, S. Prelovsek, and P. Singer, Eur. Phys. J. C6, 
471 (1999). 

K. Abe et at, Phys. Rev. Lett. 92, 101803 (2004), the 
published result has been rescaled using the latest values 
from [15]. 

D. M. Asner et al, Phys. Rev. D58, 092001 (1998). 
B. Aubert et ai, Nucl. lustrum. Meth. A479, 1 (2002). 
S. Agostinelh et al., Nucl. lustrum. Meth. A506, 250 

(2003) . 

Unless explicitly stated otherwise, charge conjugate re- 
actions are included throughout this paper. 
W.-M. Yao et al. (Particle Data Group), J. Phys. G33, 
1 (2006), and 2007 partial update for the 2008 edition. 
B. Aubert et ai, Phys. Rev. D76, 011102 (2007). 
M. J. Oregha, Ph.D Thesis, SLAC-236 (1980), J. E. 
Gaiser, Ph.D. Thesis, SLAC-255 (1982), T. Skwarnicki, 
Ph.D Thesis, DESY F31-86-02 (1986). 



